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Abstract. We calculate multicomponent radiatively driven stellar wind mod- 
els suitable for central stars of planetary nebulae. Some of these stellar winds 
may be adequately modelled using one-component models, however for some 
of them multicomponent models are necessary. We obtain a range of stellar 
parameters for different types of mass-loss. 



1. Introduction 

Hot star winds are accelerated mainly by absorption of radiation in resonance 
lines. This process can be divided into two steps: 

1. transfer of light momentum to C, N, O, Fe, etc., by absorption of radiation, 
and by Thomson scattering to free electrons, 

2. transfer of obtained momentum to predominant wind component (H, He). 

Since the acceleration of different wind species may be different, their wind ve- 
locities may differ and so hot star winds have a multicomponent nature. For 
high density winds (e.g. of galactic O stars) this multicomponent nature does 
not influence the wind structure, however for low density winds multicomponent 
effects occur, for example frictional heating, decoupling of wind components, 
etc. (see Krticka & Kubat 2001). To estimate the importance of multicom- 
ponent effects for the winds of central stars of planetary nebulae we calculate 
multicomponent wind models for these stars. 



2. Model assumptions 

Basic assumptions of our models are the following: 

• we assume a stationary spherically-symmetric flow, 

• we solve the continuity, momentum and energy equations for each compo- 
nent of the flow, namely for absorbing ions, nonabsorbing ions (hydrogen 
and helium), and electrons (see Krticka & Kubat 2001), 

• we assume solar chemical composition, 

• line radiative force is calculated in the CAK approximation (Castor, Ab- 
bott & Klein 1975), we neglect wind instabilities and magnetic fields. 
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3. Multicomponent model equations 

For each component a of the flow (i.e. accelerated ions, passive component (hy- 
drogen and helium) and electrons) we solve continuity, momentum and energy 
equations in the form of 

(r 2 PaVro) = 0, (1) 
dv ra rad 1 d , 2 \ . Qa „ V - ^ fric / \ 

«r.-dT=f- -9-- Tr M + -E + ^ 9ab , (2) 
3 p a dT a a\p a d . 2 



-rkv ra - - 

2 m n dr 



+ ^ = + E (03 + Qat) . (3) 



where p a is the density of a component a, ty a is radial velocity, a a is the 
isothermal sound speed, g a ad is the radiative acceleration either due to the line- 
transitions or due to the electron scattering, E is the electric polarisation field, 
g a r b lc is the frictional acceleration, T a is the temperature, Q T ^ d is the radiative 
heating/cooling term (calculated using electron thermal balance method, Kubat 
et al. 1999), Q c a \ is the heat exchange and Q l * b lc is the frictional heating. 

4. The frictional acceleration 

The frictional acceleration g a T b lc depends on the velocity difference via the Chan- 
drasekhar function (<? a " c ~ G(x ab )), 

G(x ab ) = d ^ , (4) 

ZX ab 

where <3?(x ao ) is the error-function and the relative velocity difference between 
wind components a and b is 

_ \Vrb ~ Vra\ 

•Lab — • \p) 

OLab 

For small relative velocity differences, x ab % 1, the Chandrasekhar function 
G{x ab ) is increasing, the flow is stable in this case, however for larger velocity 
differences, x ab > 1, G(x ab ) is decreasing (see Fig.l). The latter behaviour 
enables decoupling of wind components, the flow is unstable for larger velocity 
differences (Owocki & Puis 2002, Krticka & Kubat 2002). 



5. Examples of calculated wind models 

Frictional heating is important if the velocity difference is comparable with av- 
eraged thermal speed, x ab < 1. An example of frictionally heated wind model is 
given in Fig. 2. 

If the relative velocity difference between wind components is higher than 
the averaged thermal speed, x ab > 1, the wind components may decouple. An 
example of model with hydrogen decoupling is given in Fig. 3. 
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Figure 1. The run of Chandrasekhar function G(x a b). The frictional accel- 
eration is proportional to the Chandrasekhar function, g^ c ~ G(x a b)- 





Figure 2. The stellar wind model of star with parameters T c g — 100 000 K, 
M = 0.6 M© and log g = 6.46 (CGS). Velocities of wind components are 
nearly equal in this case and the stellar wind is heated in the central parts of 
the model due to the frictional heating 




Figure 3. The stellar wind model of a white dwarf with parameters T e g = 
100 000 K, M = 0.6 M© and log.g = 7.57 (CGS). The absorbing component is 
not able to accelerate hydrogen any more, hydrogen may fall back onto the 
stellar surface or may create the circumstellar shells (Porter & Skouza 1999). 
The stellar wind is significantly frictionally heated in this case. 
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Figure 4. Regions in T e g/ log g diagram with different kinds of stellar wind 
6. Regions in T e ff/log<7 diagram 

Using our wind models, for stars with M = 0.6 M Q we have derived regions in 
the T e ff/log<7 diagram with different types of stellar winds (see Fig. 4). Cor- 
responding evolutionary track by Blocker (1995) of a post-AGB stage is also 
plotted in this figure. Apparently, in the course of the post-AGB evolution also 
the stellar wind evolves. At the initial stages the stellar wind can be regarded 
as a one-component. During the subsequent stellar cooling the stellar wind be- 
comes frictionally heated and subsequently either hydrogen falls back onto the 
stellar surface or pure metallic wind exists. The coolest stars do not have any 
wind. 
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